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Abstract 

Exclusive charmless B decays to rj' and rj are studied in the standard effective Hamiltonian 
approach. The effective coefficients that take into account nonfactorizable effects in external and 
internal W^-emission amplitudes are fixed by experiment. Factorization is applied to the matrix 
elements of QCD, electroweak penguin operators and next-to-leading corrections to penguin Wilson 
coefficients, which are renormalization-scheme dependent, are included, but the resultant amplitude 
is renormalization-scheme and -scale independent. We show that the unexpectedly large branching 
ratio of — > rj'K^ measured by CLEO can be explained by the Cabibbo-allowed internal W- 
emission process b — > ccs followed by a conversion of the cc pair into the r]' via gluon exchanges 
characterized by a decay constant of order —50 MeV. Implications are discussed. The mechanism of 
the transition cc — > r/', if correct, will enhance the decay rates of i? ^ rj'K* by more than an order 
of magnitude. However, it plays only a minor role to the rare decays B r]' {ri)ii{p). The predicted 
branching ratio of B^ — > r/vr^ for positive p is marginally ruled out by experiment, indicating that 
the favored Wolfenstein parameter p is negative. 
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1. The inclusive and exclusive rare B decays to r]' have recently received a great deal of 
attention. The CLEO collaboration has reported very large branching ratios for inclusive rj' 
production 

B{B^ 7]'Xs; 2.0 GeV < Pr,' < 2.7 GeV) = (6.2 ± 1.6 ± 1.3) x 10-^ (1) 
and for the exclusive decay B rj'K |^ 

B{B^ 7]'K^) = (j.ltli ± 0.9) X 10-^ 

B(B^^iK^) = (5.3^|«±1.2) X lO-l (2) 

These results are abnormaly large; for example, the branching ratio of B^ — ^ rj'K^ is 
theoretically estimated to be of order (1 — 2) x 10~^ in the standard approach, which is 
too small by 2.5(T compared to experiment. It is natural to speculate that the unexpected 
observation has something to do with the unique and special feature of the rj' meson. 

Several mechanisms have been advocated to explain the copious rj' production in the 
inclusive process B ^ t]' + X: (i) the penguin mechanism b ^ s + g* followed by g* — >■ 
rj' + g via the QCD anomaly 0, Whether or not the standard b s + g* penguin 
alone is adequate to explain the data (1) is still under debate, (ii) the Cabibbo-allowed 
process b (cc)i + s (the subscript 1 denotes color singlet), followed by the transition 
(cc)i two gluons r]'. It was argued in ^ that the decay constant fjy^\ to be defined 
below, should be as large as 140 MeV in order to explain the inclusive rj' data solely by this 
mechanism, and (iii) the process b — *• {cc)s + s, analogous to (ii) but with the cc pair in a 
color-octet state, followed by {cc)s ^ rj' + X 0. 

In this Letter we will study the exclusive charmless B decays to rj' and t] and show that the 
CLEO measurement of B^ r]'K^ is indeed larger than what expected from the standard 
approach based on the effective Hamiltonian and factorization. We then proceed to suggest 
that two of the aforementioned mechanisms, though important for inclusive t]' production, 
do not play any significant role to the two-body decay B —>■ rj'K and that a reasonable 
amount of the charm content of the t]' arising from the transition cc ^ rj', corresponding to 
If^f^^l ~ 50 MeV, suffices to account for the data. Many interesting implications due to this 
mechanism are discussed. 



2. The hadronic rare B decays and their CP violation have been studied in detail in 
^ (see also f^, 10, [ill). "^^^ relevant effective AB = 1 weak Hamiltonian is 
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with (gig2)v±A = 9i7m(1 ±75)92- In Eqs. (3) and (4), O3 
and O7 — Oio are electroweak penguin operators. 



Oe are QCD penguin operators 



The Wilson coefficients Cj(/i) in Eq. (1) have been evaluated at the renormalization scale 
/i ~ nih to the next-to-leading order. Beyond the leading logarithmic approximation, they 
depend on the choice of the renormalization scheme (for a review, see [0). The mesonic 
matrix elements are customarily evaluated under the factorization hypothesis. In order to 
ensure the and renormalization scheme independence for the physical amplitude, the ma- 
trix elements have to be computed in the same renormalization scheme and renormalized at 
the same scale as Cj(/i). In the naive factorization approach, the relevant Wilson coefficients 
for external ly-emission (or so-called "class-I" ) and internal H^-emission (or "class-II" ) non- 
penguin amplitudes are given by ai = ci + c^jl^c and 02 = C2 + Ci/Nc, respectively, with Nc 
the number of colors. However, it is well known that the factorization approximation does 
not work for color-suppressed class-II decay modes (for a review, see This means that it 

is mandatory to take into account the nonfactorizable effects, especially for class-II modes. 
For B —>■ PP or VP decays (P: pseudoscalar meson, V: vector meson), nonfactorizable 
contributions can be lumped into the effective parameters ai and 02 [0, p^: 



,eff 



Ci(/i) +C2(/i) + Xi(/^) 



,efr 



C2{fi) + CM ( — +X2(/^) 
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(5) 



where Xi are nonfactorizable terms and receive main contributions from the color-octet cur- 
rent operators. Without loss of generality, the effective parameters a^^ include all the contri- 
butions to the matrix elements and the nonfactorizable terms Xiif^) take care of the correct 
/i dependence of the matrix elements fl^. Schematically, we have 



(c(/i)0(/i))^a'=f^(0*'-^^)fact, 



(6) 



where (0*'"'^*^)fact is the matrix element of the tree operator O^'^'^'^ evaluated using the factoriza- 
tion method. Therefore, a^^ are renormalization-scheme and -scale independent. Although 
nonfactorizable contributions are nonperturbative in nature and thus difficult to estimate, 
the effective coefficients af^ can be determined from experiment. From B D^*^Tr(p) and 
B J/tpK^*^ decays, we found that ai ~ (0.96 - 1.10) ^ and 02 = 0.26 ~ 0.30 P]. In 
the present analysis we assign the values 



af 



1.02, 



af 



0.28, 



(7) 



for exclusive two-body decays of B mesons. By comparing (5) with (7), we learn that 
factorization works reasonably well for class-I decay modes, but nonfactorizable effects play 
an essential role for class-II channels. 



3 



Unlike the effective coefficients a^;^, which can be empirically determined from experi- 
ment, a priori it is not clear how reliable is the factorization approach for the matrix elements 
of the QCD and electroweak penguin operators. Recently CLEO II has reported preliminary 
results of charmless hadronic B decays, namely, B —>■ K'^tt^, B — > K^n'^ 19|, dominated 
by the penguin mechanism. It turns out that the predictions based on factorization are in fair 



agreement with data Although at present we are not able to determine the analogous 
effective penguin coefficients that include nonfactorizable effects from experiment, we can 
nevertheless construct penguin Wilson coefficients that are independent of the choice of the 
renormalization scale and the renormalization scheme (for a review, see [0]). The one- loop 
QCD and electroweak corrections to matrix elements are parametrized by the matrices rh^ 
and rhe, respectively [|^, 



As a consequence. 



I + — msifi) + — me /i 



(9) 



One-loop penguin matrix elements of the operators Oi^2 have been calculated with the results 
p2| , p^] (penguin diagrams for the matrix elements of penguin operators are calculated in 

my- 
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-K + G{mc, k, /i) 



(10) 



where k parametrizes the renormalization scheme dependence. For example, for the naive 
dimension regularization (NDR) scheme and the 't Hooft-Veltman (HV) scheme, k is given 

by m 



-1 (NDR), 
(HV). 



(11) 



The function G{mc, k,fi) in Eq. (10) defined by 



G{mc, k, jj) = —A I dx x{l — x) In 
Jo 



ml — k'^x{l — x) 



(12) 



has the analytic expression: 
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for k"^ > 4m^, where k"^ is the momentum squared of the virtual gluon. Note that the 
imaginary part of Glm, k, fi), which furnishes the desired dynamic strong phase necessary 
for generating CP asymmetries in charged B decays, arises only in the time-like penguin 
diagram with fc^ > 4m^. Eqs. (9) and (10) lead to the renormalization scale and scheme 
independent penguin Wilson coefficients: 



where 



C3 = C3(^) + ^Ps(^), 

ce = C6(/i) - 



C4 = C4(yU) - -Ps(/U), 
C7(9) = C7{9)(/i) - PeifJ'), 



C5 = C5(/i) + 

C8(10) = C8(10)5 



(14) 



57r 
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a 
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Ci(/i), 
Ci(/i) + 3C2(/i)]. 



(15) 



As noted before, the Wilson coefficients Cj(/i) have been evaluated to the next-to-leading 



order in the NDR, HV and RI (regularization independent) renormalization schemes [0, ^ 
For the purpose of illustration, we use the AB = 1 Wilson coefficients obtained in HV and 
NDR schemes at /i = 4.4 GeV, = 225 MeV and mt = 170 GeV in Table 22 of [0]. 
From Eqs. (11-15) we find that £3 — cio obtained from Cj(/i) in the HV and NDR schemes 
are numerically very similar, implying the renormalization scheme independence of Cj. The 
values of q at /c^ = m^/2 are 



C3 = 0.0200 + i0.0048, 
C5 = 0.0155 + i0.0048, 
dj = -(0.0757 + i0.0558)a, 
Cg = -(1.3648 + z0.0558)a. 



Q = -0.0515 • 
cq = -0.0565 • 
cs = 0.057 a, 
Cio = 0.264 a. 



i0.0143, 
i0.0143, 



(16) 



One can also obtain q by first evaluating the renormalization scheme independent but fi 
dependent Wilson coefficients Ci{^), and then connecting them to Cj via the relation 



I + -^(^*(/^) - ^«(/^)) + ^(^e(/^) - re(/x)) 



(17) 



where the expressions of the matrices rs^e{lA are given in Using the values of Cj(/i) 

obtained at /x = m?, = 5 GeV in , the reader can check that £3 — cio obtained in this way 
are numerically in accordance with (16) except for a small discrepancy with 67. This shows 
the renormalization scale /z independence of the Wilson coefficients q. 

Two remarks are in order, (i) One can equally use the Wilson coefficients Cj(/i) or the 
renormalization-scheme independent Q(/i) to calculate the physical amplitude. However, 
the matrix elements of four-quark operators must be renormalized at the same scale /i and 
evaluated in the same renormalization scheme for the latter. The great advantage of working 
with the renormalization-scale and -scheme independent and C3— Cio is two-fold: First, the 
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corresponding tree-level operators are scale and scheme independent and hence they can be 
evaluated using factorization. Second, long-distance nonfactorizable effects for external and 
internal VT-emissions are already included in the effective coefficients a^*|. (ii) It is instructive 
to compare Cj with the leading-order Wilson coefficients Cj(/i) evaluated at /i = 4.4 GeV and 



Ag) = 225 MeV (see Table 22 of 01) 



C3 = 0.014, C4 = -0.030, C5 = 0.009, ce = -0.038, 

C7 = 0.045 a, C8 = 0.048 a, 09 = -1.280 a, Cio = 0.328 a. (18) 

It is evident that next-to-leading order corrections to the Wilson coefficients of QCD penguin 
operators are quite sizable. 

3. We are ready to employ the weak Hamiltonian (1) and the factorization approximation 
to calculate exclusive charmless rare B decays. We take the decay — > ri'K~ as an example. 
Factorization leads to 



[■(cc) / r{uu)\ 



04 + aio + 2(a6 + as) 
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X; 



X2S 
(19) 



where = C2i + -^C2i-i, 021-1 = C2i-i + ■^C2i {Nc = 3) for i >2, Xi are factorizable terms: 



X, ^ (ir-|(s«),_,|0)(V|(w6),_Ji?-)=^/^(m^-m^,)i^o'"''«), 

X2, = (Vlfe)v-JO)(i^-|(.-6),_Ji?-)=./f)(m|-m^)Fo^^(mJ,), 
X3 = WK-\{su),_jO)mub),,jB-), 



(20) 



and Ff) is the form factor defined in p6|. In the derivation of Eq. (19) we have applied the 
isospin relation X2d = X2U and equations of motion for {S — P){S + P) penguin matrix 
elements. The wave functions of the physical 77' and rj states are related to that of the SU(3) 
singlet state tiq and octet state % by 



7]' = riQ sin + rjo cos 9, V = Vs cos 9 — rjo sin 9 



(21) 



with 9 ^ —20°. When the t] — rj' mixing angle is —19.5°, the t]' and rj wave functions have 
simple expressions M 



\r]') = — + dd + 2ss), \r]) = --j=\uu + dd — ss). 



(22) 
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For simplicity we will employ 9 = —19.5° in ensuing discussion. In SU(3) limit, the t]' and 
r] decay constants become [0] 

jiuu) ^ fidd) ^ _fi,s) = fj^= 77 MeV. (23) 



The measured physical decay constants and are fairly close to /^r [0. As a result, 
in SU(3) limit, = 2X2^ for rj' production and = —X2U for the 77. In the limit of 
SU(3)-flavor symmetry, we also have 

Fo^^(O) = v^Fo^^'(O) = Ft\0), (24) 

where the factor of ^/6 comes from the normalization constant of the i]' wave function [see 
Eq. (22)]. 

For non-penguin external and internal H^-emission contributions to rj'K^, we 

employ the effective coefficients a^^ [Eq. (7)] since they take into account nonfactorizable ef- 
fects. In Eq. (19) X3 stands for the ly-annihilation contribution. In the penguin mechanism, 
X3 arises from the space-like penguin diagram. It is common to argue that 14^-annihilation 
is negligible due to helicity suppression, corresponding to form factor suppression at large 
momentum transfer, = (for a recent study, see [^). However, we see from Eq. (19) 
that it is largely enhanced in the space-like penguin diagram by a factor of m\/ {mi^nis). 
Unfortunately, we do not have a reliable method for estimating VT-annihilation, though re- 
cent PQCD calculations suggest that space-like penguins are small compared to the time-like 
ones ra. 



In terms of the Wolfenstein parametrization [^, the relevant quark mixing matrix ele- 
ments for B rj'K are 

V^,V:, = AX\p-tr^), V,,v:, = AX\l-^), Vt,V* = -AX', (25) 

where A = 0.22 and A = 0.804 for \Vcb\ = 0.039. Since the parameters p and rj have not been 
measured separately, it is desirable to have the results presented for allowed regions of p and 
rj. As p and r] are approximately constrained by 0.2 < r] < 0.4 and —0.3 < p < 0.3 (see, e.g., 
|12[), we will fix ~ 0.30 and take p = 0.30 as well as p = —0.30. We also need input of 
other parameters. For quark masses we use m„ = 5 MeV, rrid = 10 MeV, = 175 MeV, 
rric = 1.5 GeV and rrib = 5.0 GeV. Heavy-to-light mesonic form factors have been evaluated 
in the relativistic quark model [^ To be specific, we use F^^ (0) = 0.254/v^ [^ Qand 
F(f^(0) = 0.34 For the lifetime of B mesons, the world average is r(Erf) = (1.55 ±0.04) 
ps and t{B^) = (1.66 ± 0.04) ps |Q. Since what CLEO has measured is the combined 
branching ratio 

B{B^ r]'K^) = i [B{B+ r]'K+) + B{B- , (26) 



^It should be stressed that the form factors given in ||22l for B{D) vr*', rj, rj' transitions do not take 
into account the normaUzation constant of the neutral meson wave functions. 



7 



we shall present the results for the combined one. Neglecting the l^-annihilation contribution 
characterized by the parameter X3, we find that in the absence of the transition cc — > rj' (see 
Table I) 



e(S± ^ ^'i^±) = (1.4 _ 1.8) X 10-^ for 4f ^ = 0, (27) 

which is substantially smaller than the CLEO measurement (2). It is worth emphasizing 
that if the Wilson coefficients are evaluated only to the leading order [cf. Eq. (18)], then 
the resultant predictions become even worse: B{B^ ^ v'K^)l.o. = (0.5 - 0.8) x 10"^ In 
fact, the leading order calculation is already ruled out by the recent CLEO measurement of 
n^K and n^K^ pO]. 



We briefiy discuss the uncertainties associated with the prediction of B{B^ rj'K^) and 
the numerical results given in Table I. We have neglected ly-annihilation, space-like penguin 
diagrams and nonfactorizable contributions to the matrix elements of penguin operators; all 
of them are difficult to estimate. Other major sources of uncertainties come from the choice 
of form factors and light quark masses. It is natural to ask if the CLEO result for B^ rj'K^ 
can be explained by different choices of the strange quark mass rris and form factors. For a 
given set of form factors, it is easily seen that the branching ratio increases with decreasing 
m<j. For example, we find a large branching ratio i3(i?^ — >■ rj'K^) = 6.5 x 10~^ for = 60 
MeV and p < 0. However, the use of a small strange quark mass is in confiict with a recent 
measurement in r decay, msirriT-) = 212II5 MeV 0]. It is also true that the measured 
branching ratio of the r]'K mode can be accommodated by choosing large form factors, e.g., 
F^^{0) ^ VGFo'^'iO) ^ 0.61. However, this will break the SU(3)-symmetry relation (24) 
very badly. All the existing QCD sum rule, lattice and quark model calculations indicate 
that F(f^ (0) ^ 0.33 or less (for a review, see [p5]). The recent observed charmless decays 
B — > K'^n^, B K^vr^ dominated by the penguin diagram are governed by the form 
factor Ff^'"^ . We find that the predictions for B Kn based on the standard approach 
with F(f^ (0) = 0.33 are in agreement with data. Therefore, in view of the SU(3) relation 
(24) for form factors, it is very unlikely that F^^{0) and 76^0^''' (0) can deviate much from 
0.33. We should also stress that it is dangerous to adjust the form factors in order to fit a 
few particular modes; the comparison between theory and experiment should be done using 
the same set of form factors for all channels. We thus conclude that the large branching 
ratio of B^ —>■ rj'K^ observed by CLEO cannot be explained by the conventional approach. 

The sizable discrepancy between theory (27) and experiment (2) strongly suggests a new 
mechanism unique to the t]' meson. It appears that the most natural one is the Cabibbo- 
allowed process b — > ccs followed by a conversion of the cc pair into the t]' via two gluon 
exchanges [see Eq. (19)]. This new internal PF-emission contribution is important since 
its mixing angle VcbV*^ is as large as that of the penguin amplitude and yet its Wilson 
coefficient a^^ is larger than that of penguin operators. The decay constant f^f^\ defined 

as (0|c7^75c|?7') = if^f^\i_i, can be reduced via the OPE to a matrix element of a particular 

^The prediction B{B^ -q' K^) = 3.6 x 10^*^ given in Q is too large by about a factor of 2 because the 
normalization constant (i.e. 1/V6) of the 77' wave function was not taken into account in the form factor 
Fq^ (see also the previous footnote). 
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dimension-6 pseudoscalar gluonic operator It has been estimated in |]2S|, based on the 
OPE, large- iVc approach and QCD low energy theorems, that |/?^^| = (50-180) MeV. (An 



independent estimate of f^f^^ based on the instanton-liquid model is given in |3^.) It was 



claimed in p8| , |^ that If^y'^ \ ~ 140 MeV is needed in order to exhaust the CLEO observation 
of r]'K^ and B rj' + X hj the mechanism b^cc+s^rj' + s via gluon exchanges. 

In view of the fact that /^r""* is only of order 50 MeV, a large value of f^?'^^ seems to be 
very unlikely. We believe that it is more plausible to take the lower side of the theoretical 
estimate, namely l/^f'^^l ~ 50 MeV. On the other hand, it is claimed in that If^f^'^l < 40 
MeV. From Eq. (19) it is clear that this additional contribution will enhance the decay rate 
of B^ — » rj'K'^ provided that the sign of f^?^^ is opposite to that of /^r""*. (Note that the 
estimate of /^^^ in does not fix its sign.) For /^f^-* = —50 MeV, we obtain 

B{B^^r]'K^) = (5.8 -6.7) X 10-^ 

S(5° ^ v' = (5.7 - 5.9) X 10-^ (28) 

which are in agreement with the CLEO measurement (2), especially when p is negative. 



Contrary to what has been claimed in |28], it is not necessary to introduce a large value of 



f^f'^'' to account for the experimental observation of B^ rj'K^. 

Suppose the dominant mechanism for the observed large inclusive B^ — * 77' + X signal 
comes either from the b ^ s g* penguin followed by the transition g* gt]' via the QCD 
anomaly, or from the process b — > {cc)s + s, followed by the conversion {ccjs — > rj' + X. 
Will these mechanisms still play an important role in exclusive B^ rj'K^ decays ? It 
is easily seen that since the above two mechanisms involve a production of a gluon before 
hadronization, they will not make contributions to two-body exclusive decays unless the 
gluon is soft and absorbed in the wave function of the 77'. Another possibility is the penguin- 
like process b —>■ s + g*g* —>■ s + rj', but it is a higher order penguin mechanism. 

The factorizable contribution to B^ — > rjK' has a similar expression as (19) for B^ — >■ 
ri'K~ with the replacement 

/f)-/f)^-sin^/f\ (29) 

where 6 ~ —20° is the rj — rj' mixing angle. Because of the destructive interference in the 
penguin diagrams due to the relation X2S = —X2U in SU(3) limit and because of the smallness 
of the cc rjo contribution to rj production, the decay B rjK is suppressed. We see from 
Table I that the branching ratio of B ^ rjK is smaller than that of B ^ rj'K by an order of 
magnitude. Nevertheless, a measurement of B{B — »• rjK) of order (2 — 5) x 10~^ will confirm 
the importance of the cc — > rj' mechanism. From Table I we see that the electroweak penguin 
effects are in general very small, but they become important for B rjK and B rjK* 
decays due to a large cancellation of QCD penguin contributions in these decay modes. It is 
worth stressing that CP asymmetries in B"^ r]{ri')K^{K*^) decays will be largely diluted 
in the presence of the cc conversion into the rj' p6|. 



•^The mixing angle of the tj' with the cc state is of order 7° for ~ 50 MeV. 
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The mechanism of cc conversion into 77', if correct, is probably most dramatic in the decay 
B rj'K* . The factorization amphtude for ri'K*~ is given by 



(afX[ + af XL + afX'^) + V^,V:,af X'.^ff^ ' ^^^^^ 



/f' 



[04 + aio)X[ + ( 2a3 - 205 - ^07 + ^ag ) X!^^ 



ml, 



+ 04 + aio — 2(a6 + cis) 



m 



B 



{rus + mu){mb + rUu) 



X' 



(30) 



where 



{K*-\{-su),_^Q){i\{uh),_,\B-) = 
(r/'|(gg),_JO)(ir*-|(.-6),_Ji?-) = 
(Vir*-|(sn),_jO)(0|(n6),_Ji?-). 



-2fK'mK'F^'^\m\,){e-pB), 



(31) 



Note that, contrary to i? ^ rj'K decay, there is no contribution proportional to X[ arising 
from the {S — P){S + P) part of the penguin operators and Os- The dependence of the 
form factors Fq, Fi, Aq, defined in |^^, can be calculated in the framework of the relativistic 
quark model. A direct calculation oiB ^ P and B form factors at time-like momentum 
transfer in the relativistic light-front quark model just became available recently [0. It is 
found that Aq, Fi exhibit a dipole behavior, while Fq shows a monopole dependence. At 



g = 0, the relevant form factors are |]32 | 



Ff^(O) = 0.34, Ff"(0) = 0.26, A^^* (0) = 0.32, A^^O) = 0.28. 



(32) 



We then utilize the above-mentioned behavior to compute the form factors at the desired 
momentum transfer squared. Owing to the small penguin contributions, the standard ap- 
proach's estimate of B{B rj'K*) is small, of order (2 — 5) x 10~^ for B^ rj'K*^ and 
1 X 10~^ for B^ —>■ rj'K*^ . Since the branching ratio due to the process 6— >cc + s— ^r^' + s 
is of order 10^^, B{B rj'K*) is greatly enhanced to the order of (1 — 2) x 10~^. It 
has been argued in that B{B rj'K*) is about twice larger than that of -B ^ ^'K, 
which is certainly not the case in our calculation. It is interesting to note that the ra- 
tio B{B r]'K*)/B{B rjK*) is much less than unity in the naive estimate but it 
becomes greater than unity in the presence of cc conversion into the tiq (see Table I). 



For B — i> ri'{ri)7T{p) decays, the mechanism of cc r]o is less dramatic since it does not 
gain mixing-angle enhancement as in the case of -B — > r]'{ri)K{K*). We see from Table I that 
the full branching ratios in general are close to the predictions in the conventional approach, 
indicating the minor role played by the charm content of the rj'. Three remarks are in order. 
First, there is a large enhancement of order / {mqimb) {q = u, d) occurred in the matrix 

elements of {S — P){S + P) penguin operators. The calculation in this case is thus sensitive 
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Table I. Averaged branching ratios for charmless B decays to rj' and r] [see Eq. (26)], where "Tree" 
refers to branching ratios from tree diagrams only, "Tree+QCD" from tree and QCD penguin 
diagrams, and "Full" denotes full contributions from tree, QCD and electroweak (EW) penguin 
diagrams in conjunction with contributions from the process cc — > tjq. Predictions are for k'^ = 
ml/2, /^r^ = —50 MeV, rj = 0.30, p = 0.30 (the first number in parentheses) and p = —0.30 (the 
second number in parentheses). 



Decay 


Tree 




Tree+QCD 




Tree+QCD+EW 


Full 




Exp. |1,3 


B± - 




1.45 X 10" 


-7 


(1.40, 


1.82) 10" 


-5 


d 35 1 76) 10" 


-5 


(5.83, 6.67) 10" 


-5 


(7.1^2:? ±0.9) 10-^ 


B± - 




4.12 X 10" 


-7 


(2.73, 


6.54) 10" 


-7 


ro 21 1 05) 10" 


-6 


(1.94, 5.15) 10- 


-6 


< 0.8 X 10" 


-5 


B±- 


1 


2.10 X 10" 


-7 


(1.57, 


5.23) 10" 


-1 


(1.67, 4.57) 10" 


~7 


(1.14, 1.62)10" 


-5 


< 29 X 10- 


-5 


B±- 




6.23 X 10" 


-7 


(0.41, 


2.09) 10" 


-6 


(0.56, 2.79) 10" 


-6 


(2.95, 7.79) 10" 


-6 


< 24 X 10- 


-5 


B±- 


tj'tt^ 


1.88 X 10" 


-6 


(8.94, 


8.13) 10" 


-6 


(8.91, 8.09)10" 


-6 


(1.36, 1.24) 10" 


-5 


< 4.5 X 10" 


-5 






5.34 X 10" 


-6 


(9.19, 


1.89) 10" 


-6 


(9.32, 1.91)10" 


-6 


(1.07, 0.22) 10" 


-5 


< 0.8 X 10" 


-5 


i3±- 


■q'p^ 


3.90 X 10" 


-6 


(0.45, 


2.37) 10" 


-5 


(0.45, 2.39) 10" 


-5 


(0.47, 1.83) 10" 


-5 






B±- 


■qp^ 


1.16 X 10" 


-5 


(1.05, 


1.84) 10" 


-5 


(1.06, 1.80) 10" 


-5 


(1.15, 1.65) 10" 


-5 






Bd-^ 




7.93 X 10" 


-9 


(1.45, 


1.54) 10" 


-5 


(1.40, 1.49)10" 


-5 


(5.73, 5.92) 10" 


-5 


(5.31^2 ±1-2) 


10-^ 


Bd^ 




1.58 X 10" 


-8 


(0.31, 


1.07)10" 


-7 


(1.52, 3.14)10" 


-7 


(2.67, 3.29) 10" 


-6 






Bd^ 




6.05 X 10" 


-9 


(0.93, 


1.50)10" 


-7 


(0.73, 1.17)10- 


-7 


(1.28, 1.36) 10" 


-5 


< 9.9 X 10" 


-5 


Bd^ 




1.32 X 10" 


-8 


(4.56, 


6.88)10" 


-7 


(0.82, 1.13)10" 


~6 


(4.17, 4.86)10" 


-6 


< 3.3 X 10" 


-5 


Bd^ 




0.34 X 10" 


-9 


(2.27, 


6.81)10" 


-6 


(2.20, 6.61)10" 


-6 


(3.72, 9.19)10" 


-6 


< 2.2 X 10" 


-5 


Bd^ 




1.98 X 10" 


-9 


(0.73, 


2.43) 10" 


-6 


(0.71, 2.36) to- 


-6 


(0.96, 2.83) 10" 


-6 






Bd^ 


r]'p" 


1.03 X 10" 


-8 


(1.57, 


4.18) 10" 


-6 


il. 67, 4.45) 10" 


'6 


(0.77, 2.78) 10" 


-6 






Bd^ 




5.86 X 10" 


-8 


(2.84, 


2.60) 10" 


-7 


(3.20, 3.24) 10- 


-7 


(1.67, 1.74)10" 


-7 


< 8.4 X 10" 


-5 



to the light quark masses and rrid- Second, the branching ratio of rjir"^ is predicted 

to be 1 X 10^^ for positive p, while experimentally B{B^ rjh^) < 8.0 x 10""^ for h = n, K 
This again indicates that a negative p is preferable. Third, the decays rjp^ are 

overwhelmly dominated by tree diagrams. 

4. To summarize, in view of the recent unexpectedly large branching ratios for inclusive 
and exclusive charmless B decays to r/', we have analyzed the exclusive B decays to rj' and 
7] in detail. We showed that the process b^cc + s^rj' + s via gluon exchanges with 
f^f^^ ~ —50 MeV is adequate to explain the discrepancy between theory and experiment. 
A measurement of B{B rj'K*) of order 10~^ will give a strong support of this gluon 
mechanism of rj' production. Moreover, we found that the decay rate of B ^ rj'K* is larger 
than B — > rjK*, contrary to what expected from the standard approach. The predicted 
branching ratio of B^ 777?^ for positive p is marginally ruled out by experiment, implying 
that a negative p is preferable. 

ACKNOWLEDGMENT: This work was supported in part by the National Science Council 
of ROC under Contract Nos. NSC86-2112-M-001-020 and NSC86-2112-M-001-010-Y. 

Note added: The exclusive decays of B to rj' and rj were also discussed recently by A. Ali 
and C. Greub (|hep-ph/9707251|) , A. Datta, X.G. He and S. Pakvasa (|hep-ph/ 9 7072591) with 



conclusions different from ours. 
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